Pseudo-zero-mode Landau levels and collective excitations in bilayer graphene 
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Bilayer graphene in a magnetic field supports eight zero-energy Landau levels, which, as a tunable 
band gap develops, split into two nearly-degenerate quartets separated by the band gap. A close 
look is made into the properties of such an isolated quartet of pseudo-zero-mode levels at half filling 
in the presence of an in-plane electric field and the Coulomb interaction, with focus on revealing 
further controllable features in bilayer graphene. The half-filled pseudo-zero-mode levels support, 
via orbital level mixing, charge carriers with nonzero electric moment, which would lead to field- 
induced level splitting and the current-induced quantum Hall effect. It is shown that the Coulomb 
interaction enhances the effect of the in-plane field and their interplay leads to rich spectra of 
collective excitations, pseudospin waves, accessible by microwave experiments; also a duality in the 
excitation spectra is revealed. 

PACS numbers: 73.43.-f,71.10.Pm,77.22.Ch 
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I. INTRODUCTION 

Graphene, a monolayer of graphite, attracts a great 
deal of attention, both experimentally^^ and theoret- 
ically jiiSA 7 ^ for its unusual electronic transport, char- 
acteristic of "relativistic" charge carriers, massless Dirac 
fermions. Dirac fermions give rise to quantum phenom- 
ena reflecting the particle-hole picture of the vacuum 
state, such as Klein tunneling^ and, especially in a mag- 
netic field, such peculiar phenomen a 10 ' 11 ' 12 ' 13 ' 14 as spec- 
tral asymmetry and induced charges, which are rooted 
in the chiral anomaly (i.e., a quantum conflict between 
charge and chirality conservations). Graphene provides 
a special laboratory to test such consequences of quan- 
tum electrodynamics. Actually, the half-integer quantum 
Hall (QH) effect and the presence of the zero-energy Lan- 
dau level observed 1 ^ in graphene are a manifestation of 
spectral asymmetry. 

Bilayer graphene is equally exoti o 15 ' 16 ' 17 as monolayer 
graphene. It has a unique property that the band gap 
is controllabl o 18 ' 19 ' 20 ' 21 ' 22 by the use of external gates or 
chemical doping; this makes bilayers richer in electronic 
properties. In bilayer graphene intcrlayer coupling mod- 
ifies the intralayer relativistic spectra to yield quasipar- 
ticles with a parabolic energy dispersionJ£ The particle- 
hole structure still remains in the "chiral" form of a 
Schroedinger Hamiltonian, and there arise eight zero(- 
energy)-mode Landau levels (two per valley and spin) in 
a magnetic field. 

Zero-mode Landau levels, specific to graphene, deserve 
attention in their own right. They would show quite 
unusual dielectric respons e 23 ' 2 while they carry normal 
Hall conductance e 2 /h per level; for bilayer graphene di- 
rect calculations 25 indicate that the zero-modes show no 
dielectric response for a zero band gap but the response 
grows linearly with the band gap. In bilayer graphene, 
with a tunable band gap, the zero-mode levels split into 
two quartets separated by the band gap at different val- 
leys. Such an isolated quartet of "pseudo" -zero- mode 
levels remains nearly degenerate and, as noted earlier j^ 5 - 



the level splitting is enhanced or controlled by an in- 
plane electric field or by an injected current; this opens 
up the possibility of the current-induced QH effect for 
the pseudo- zero-mode sector around half filling, i.e., at 
filling factor v — ±2. 

The purpose of this paper is to further examine the 
properties of the pseudo-zero-mode levels, especially, co- 
herence and collective excitations in the presence of an 
external field and the Coulomb interaction. The pseudo- 
zero-mode levels at half filling support, via orbital level 
mixing, charge carriers with nonzero electric dipole mo- 
ment, which is responsible for field-induced level splitting 
and the current-induced QH effect. Along this line our 
discussion comes in contact with the works of Barlas et 
al2& and Abergel et al^L who, from the viewpoint of QH 
ferromagnets, studied the interaction-driven QH effect in 
the nearly-degenerate octet of zero-mode levels in bilayer 
graphene. Our paper partly extends their analysis by re- 
vealing an interesting interplay between the in-plane field 
and Coulomb exchange interaction, which leads to rich 
spectra of collective excitations, accessible by microwave 
experiments. We shall find a duality in the excitation 
spectra and, under certain circumstances, an instability 
in pseudospin textures. 

In Sec. II we briefly review some basic features of the 
pseudo-zero-mode levels in bilayer graphene. In Sec. Ill 
we construct a low-energy effective theory for the half- 
filled pseudo-zero-mode sector. In Sec. IV we examine 
the spectrum and collective excitations in it. In Sec. V 
we study the microwave response of collective excitations. 
Section VI is devoted to a summary and discussion. 



II. BILAYER GRAPHENE 

Bilayer graphene consists of two coupled hexagonal lat- 
tices of carbon atoms, arranged in Bernal A'B stack- 
ing. The electron fields in it are described by four- 
component spinors on the four inequivalent sites (A, B) 
and (A',B') in the bottom and top layers, and their 
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low-energy features are governed by the two inequivalent 
Fermi points K and K' in the Brillouin zone. The in- 
tralayer coupling 70 = jab is related to the Fermi veloc- 
ity v = (V3/2) a^o/h w 10 6 m/s (with a L = 0.246nm) 
in monolayer graphcnc. The interlayer couplings 71 = 
jA' b and 73 = jab' are one-order of magnitude weaker 
than 70; numerically, 28 71 w 0.30 eV, 73 « 0.10 eV and 
70 ~ 2.9 eV. 

Actually, interlayer hopping via the (A' , B) dimer sites 
modifies the intralayer relativistic spectra to yield spec- 
tra with a quadratic dispersion and the characteristic 
cyclotron energy lu c = 2v 2 /(ji£ 2 ) ~ 3.9 7?[T] meV, 
with a magnetic field B[T] in Tesla. The low-energy 
branches, thereby, are essentially described by two- 
component spinors on the {A, B') sites (with the high 
energy branches being separated by a large gap ~ 71). 
The effective Hamiltonian is written as^ 



H 



d 2 x 



^(H+- eA )i> + x t (H--eA ) X 



H^ = Ho + Has, 



Ho = u c 



-a 2 + A at 



-(at) 2 + Aa 



H a 



U fl-zah 



-(1 — z cur) 



(2.1) 



together with coupling to electromagnetic potentials 
(Ai,Ao). Here, assuming placing graphcnc in a uni- 
form magnetic field B > 0, we have rescaled the ki- 
netic momenta IF = — idi + eAj with the magnetic 
length 1 = 1/ 

and aJ = \/2eB (H x + iFI y ) so that [a, o^] = 1 
Ai — > B (— y, 0) to supply a strong magnetic field B nor- 
mal to the sample plane. The field ip — (ipA, V's 
to the K valley with H+ — H^=i while \ — (xb 



V2^B (H x - o^y 



/IF,) 

we set 
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refers to the K' valley with H- = H^=-i. For simplic- 
ity we ignore weak Zeeman coupling and suppress the 
electron spin indices. 

In Ho the linear kinetic term, leading to "trigonal 
warping" , represents direct interlayer hopping via 73 , 
with A = £ (73/70X^2 vo/£)/u c w ±0.3/v/b[T]. 

The Has takes into account a possible layer asymme- 
try, caused by an interlayer voltage AAq, which leads 
to a tunabl o 18 i 19 ' 20 gap U w eAAg between the con- 
duction and valence bands. The O(za'a) terms in Has 
represent a kinetic asymmetry related to the charge de- 
pleted from the dimer sites; note that it is very weak, 
with z = 2lu c /ji « 0.026 x B[T] < 1. 

While the linear spectra are lost, the bilayer Hamilto- 
nian still possesses a key feature of relativistic field the- 
ory, the particle-hole (or chiral) structure of the quantum 
vacuum: When the tiny O(zU) asymmetry is ignored, 
the spectrum of H^ , in general, is symmetric about zero 
energy e = 0, apart from possible e = ±{7/2 spectra that 
evolve from the zero-energy modes of Ho- Indeed, for 
A = the spectrum of H^ consists of an infinite tower 
of Landau levels \n,yo) of paired positive and negative 



energies, 

e n = s n w cy /\n\(\n\ - 1) + {U/2lu c ) 2 -\^zU, (2.2) 

labeled by integers n = ±2, ±3,..., and p x (or yo = 
£ Px)l s n = sgn{n} = ±1 specifies the sign of e n . 

In addition, there arise nearly-degenerate Fandau lev- 
els carrying the orbital index |n| = and \n\ = 1, with 
spectrum 



-\n\=0 



tU/2, 



\n\ = l 



C (?7/2) (I-0). (2.3) 



From now on we take, without loss of generality, U > 
and denote n — 0± and n = ±1 to specify these pseudo- 
zero-mode levels. There are four such pseudo-zero-mode 
levels (or two levels per spin) at each valley and they 
reside on different layers; the (0+,l) quartet on A sites 
at £ = 1 valley is separated from the the (0_, —1) quar- 
tet on B' sites by a band gap U. Each quartet remains 
degenerate, apart from 0(z U) fine splitting. 

The presence of the pseudo-zero-mode levels and their 
double-fold degeneracy (per spin and valley) both have a 
topological origin, and are consequences of spectral asym- 
metry, or the nonzero index of the Hamiltonian H^\u~*0: 



Index[ft f |(7^o] = J d 2 x2 x (eB/2n), 



(2.4) 



which is tied to the chiral anomaly in 1+1 dimensions. 
This degeneracy is unaffected by the presence of trigonal 
warping A ^ alone^ but is affected by nontrivial diag- 
onal components in H^. Indeed, the kinetic asymmetry 
~ zo^a leads to tiny level splitting and an electric field 
~ diAo can also enhance the splitting. In other words, 
the pseudo-zero-mode levels have an intrinsic tendency 
to be degenerate, but this at the same time implies that 
their fine structure or the way they get mixed depends 
sensitively on the environment. 

The main purpose of the present paper is to study such 
controllable features of the isolated pseudo-zero-mode 
quartet in the presence of external fields and Coulomb 
interactions. We are particularly interested in the prop- 
erties of such a quartet at half filling, where mixing of 
the zero-modes leads to nontrivial coherence effects and 
collective excitations. For definiteness we focus on the 
n = (0+, 1) sector at £ = 1 valley, i.e., around filling fac- 
tor v = 2 (or v — 1 when the spin is resolved) , and ignore 
the presence of other levels which are separated by rela- 
tively large gaps; the v = —2 case is treated likewise. We 
ignore the effect of trigonal warping, which causes only 
a negligibly small level splitting 2 ^, of 0(z\ A ) < 1/10 3 , 
apart from a common level shift of 0(z\ 2 ). 

To project out the n = (0+,l) sector let us make 
the Fandau level structure explicit by the expansion 
i>(x,t) = *Z, n ,y (x\ii,yo) il>n(yo,t), with the field opera- 
tors obeying {ip m (yo, t),ipl(y ,t)} = S mn S(y - y' ). The 
charge density p- p (t) — J d 2 xe !p ' x ifj^ip is thereby writ- 
ten as 



P-p = 7p 

k.n— — oo 



(2.5) 
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where 7 P = e~ rp2 / 4 ; r = (r x ,r y ) — (i£ 2 d/dyo,yo) stands 
for the center coordinate with uncertainty [r x ,r y ] = it 2 . 
For the (0+,l) sector the relevant coefficients are given 
by^ <?oo(p) = 1, 9ii (P) = l-*V/2, ffio(p) = Up/y/2, 
and goi (p) = itpW^, with p = p y + ip x . 

Let us put the ipo + (yo,t) and ipx(yo,t) modes into a 
two-component spinor ^> = (ipo + ,tpi) t , and define the 
pseudospin operators in the (0+, 1) orbital space, 

S* p = 7 P / dy &±v»e i »-'% (2.6) 

where <r p = (l,cr a ) (/x = ~ 3) with Pauli matrices a a 
and <j° = 1. The charge density p- p projected to the 
(0+, 1) sector is thereby written as 

p. p = 2 (w° p S _ p + w 3 p S 3 _ p + w p Si p + w 2 p S 2 _ p ), (2.7) 

where tu° = 1 - ^ 2 p 2 /4, = £ 2 p 2 /4, w p = iipy/V^, 
and w 2 (p) — ilp x j\/2. Note that (S^_ p , ^ip)* acts as a 
vector under rotations about the z axis || J5. 

The Hamiltonian H of Eq. (|2.ip is projected into the 
(0 + , 1) sector to yield 



AH 



d 2 xeA p + dy m^a 3 



(2.8) 



where m = zU/4; the common energy shift (1 — z/2)U/2 
of the (0+, 1) sector has been isolated from AH. One 
can equally write AH as 



AH = 2J2K S -p> 
p 

with V!^ being the Fourier transform of T-^x, t), 

V° = -m-eAo - ±e£ 2 V-E, 
V 3 = m- ±e£ 2 V-E, 



(2.9) 



(V 2 ,V') = —j= (E x ,E y ), 



(2.10) 



where E|| = (E Xl E y ) = — d^A^ denotes the in-plane elec- 
tric field. 

The Coulomb interaction also simplifies via projection. 
The pseudo-zero-modes at each valley essentially lie on 
the same layer, apart from a negligibly small admixture 
of O(oj c /ji) ~ 10~ 3 . One may thus retain only the in- 
tralayer interaction for 'J, 



H l 



1^2 v p : p-pPp ^ 



(2.11) 



where v p = 27ra/(eb|p|) is the Coulomb potential with 
the fine-structure constant a = e 2 /(47reo) w 1/137 and 
the substrate dielectric constant e^; J2 P = / d 2 p/(2n) 2 . 
The normal-ordered charges in H c are rewritten as 

: P-pPp ■ = P-pPp - A, 

A = 2 1 2 p {\w£\ 2 S° + 2w° p w 3 p S 3 } (2.12) 



under a symmetric integration over p. 

The pseudospin operators Sp 1 obey the S\J(2)xW 00 al- 
gebraic 

[S;, S b k ] = c(p, k) ie abc S c p+k - i s(p, k) 5 ab 5 p+k , 

[S° p , S£] = [S;, S£] = -i a(p, fc) 5; +k , (2.13) 
where (a, 6) runs over 1 ~ 3, and 

S (,, fc )=sin(^) 



(2.14) 



p x k = e IJ pikj = p x k y — p y k x ] for c(p, k) set sin(- • •) 
cos(- • ■) in s(p, k). 



III. PSEUDOSPIN TEXTURES 

In this section we study the properties of the pseudo- 
zero-mode levels at half filling, using the projected Hamil- 
tonian H = AH+H C and the charge algebra ()2.13j) . with 
focus on orbital mixing of the zero-modes. Let us suppose 
that such a half-filled state is given by a classical config- 
uration where the pseudospin points in a fixed direction 
in pseudospin space, i.e., >!5 p= o = \ N e n a and n a n a = 1 
with the total number of electrons N e = 2 S p=0 . 

Note that n 3 = 1 corresponds to the filled n — 0+ level 
with the vacant n = 1 level while n 3 = — 1 represents the 
filled n = 1 level. The direction n = (n 1 ,n 2 ,n 3 ) would, 
in general, vary in response to the external field A , and, 
as n tilts from 



n 3 



■ I . the n = 0+ and n = 1 lev- 
els start to mix. For selfconsistency we assume that Ao 
represents a uniform in-plane electric field E|| = — d x Ao 
and that it leads to a homogeneous state \G) of uni- 
form density po — v/(2irl 2 ) (with filling factor v = 2 
for the spin-degenerate v — 2 state or v = 1 for the 
spin-resolved v = 1 state). We thus consider all classi- 
cal configurations with (G\S p=0 \G) = \ N e n a and single 
out the ground state or the associated n a by minimizing 
the energy (G\H\G). For n a we use the parametriza- 
tion n 1 = sin# cos</), n 2 = sin# sin</>, n 3 — cos 9, with 
— 7r < 6 < 7r and < <f> < tt. We denote expectation 
values (G\0\G) = (O) for short. 

Let us first substitute the charge and pseudospin, 
(S p ) = ±p 5 P A and = ±p n a 5 p , , into AH. (Here 
Sp.o = (2?r) 2 <5 2 (p), and So,o = f d 2 x equals the to- 
tal area so that N e = po^o,o-) This yields the re- 
sponse of the classical state \G) to an external probe, 
(AH) =p J d 2 yLn^V^ with 

n V-p» = _ e7 4 Q + m ( cos Q _ 1) _ E|| • d e , 

ei et 
d e = -—{ n 2 ,n l ) = --j= sine fi|| , (3.1) 

where fi|| = (n x ,h y ) = (sin cj>, cos <fi) is an in-plane unit 
vector. 
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Note that the half-filled pseudo-zero-mode state has 
an in-plane electric dipole moment d e of strength 
(e£/V2) | sin 6>| per electron, proportional to the in-plane 
component (n 2 ,-n 1 ) of the pseudospin. Mixing of the 
n = 0+ and n = 1 modes gives rise to this dipole 
and its in-plane direction fi|| is related to the relative 
phase between them; i.e., under U(l) transformations 



-i<p T ( 



ia I" 1 nil rotates, 



(S 1 + iS 2 ) -> e l<t >" 1 (S 1 + iS 2 ) or <f> -> + rc i. (3.2) 

This tells us that a change in relative phase caused by a 
change in direction of the dipole is physically observable. 

In the absence of the Coulomb interaction, n a naturally 
points in the — V a direction and, as a result, an in-plane 
field, coupled to the the electric dipole, works to enhance 
the pseudo-zero-mode level splitting, 25 



Ae 



i 2 + e 2 £ 2 ~E 2 /2. 



(3.3) 



We shall discuss below how the Coulomb interaction 
modifies this. 

The calculation of the Coulomb energy (G\H C \G) = 
(H ) requires the knowledge of pseudospin structure fac- 
tors (G\S£S»\G) = (SgSQ, which, for the present half- 
filled state with pseudospin oc n a , are given by 



1 



S>11/ 7p<Wq,o + n^n" poS p ,o t> q , 



(3.4) 



where /i and v run from to 3, with n = 1 and {n a } — 
(n 1 ,™ 2 ,™ 3 ); (Ot = with = and 



„31 



,12 



cos 2 9 cos 2 i 



sin 



cos 2 9 sin 2 



- cos <p, 
sin 9 cos 9 cos d> + i sin 9 sin < 



sin t) cos t) sin < 



i sin v cos <p, 
■f i cos 6* cos 20. (3.5) 



See Appendix A for details. Actually, the normal-ordered 
correlation functions take simpler form 

( -SgSq :) = --pon^n" (7p5 p + q ,o-po <W £q,o), (3.6) 

with which one can cast the Coulomb energy in the form 
1 



(H v ) = -^N e ^2v p e-« / 2 C e (q 2 /4), 
p 



(3.7) 



where q 2 = £ 2 p 2 and Cg(x) = 1 — x(l — x) (1 — cos6*) 2 . 
From (H c ) we have omitted a constant (N e /2) pov p ^ 
which is removed when the neutralizing positive back- 
ground is taken into account. Integrating over p yields a 
typical scale of the Coulomb exchange energy, 



a 



(3.8) 



The effective Hamiltonian H c s = (AH + H c ) is con- 
veniently written in x space as H e g = p J d 2 x7i c ff with 

n cS = -eA +E(9), 
E(9) = --Vi +m(cos9- 1) +£ sm9 

+ ^F 1 (l~cos0) 2 , (3.9) 

where £ = (el/ y/2) E|| ■ nil . In the present notation (with 
Po J <i 2 x — > N e ) H, e s stands for energy per electron in 
state \G) with pseudospin cx n. The Coulomb correla- 
tion energy oc V\ consists of a negative uniform compo- 
nent —V\/2 [relative to the zero of energy (U/2) (1 — z/2)] 
and a polarization dependent component which alone fa- 
vors 9 — 0, i.e., the filled n = 0+ level, and which varies 
continuously by an amount AE C = (1/8) V\ as n sweeps 
in pseudospin space. The Coulomb interaction thus sig- 
nificantly enhances the pseudo-zero-mode level splitting. 

The stable configuration of the half-filled zero-mode 
state \G) is determined by minimizing E(9) with respect 
to n, or (9,<f>). Obviously E(9) depends on 4> through 



£, which attains a maximum when n 



E||, or £ 



e£|En|/v2- Accordingly it is convenient, without loss 
of generality, to suppose that the in-plane field E|| and 
xi|| lie along the y axis, £ = elE y j\j2 > and = 0. 
With this choice the "1", "2" and "3" axes in pseudospin 
space coincide with the y, x and — z axes in real space, 
respectively. We adopt this choice and set (n 1 ,?^ 2 ,^ 3 ) = 
(sin 9, 0, cos 9) in what follows. 

One can now look for possible ground-state configura- 
tions by writing down a phase diagram as a function of 
m and E|| . For clarity of exposition we leave it for a later 
stage and here study collective excitations over a given 
ground state \G)\ n . We focus on a special class of low- 
energy collective excitations, pseudospin waves, that are 
rotations about the energy minimum \G)\ n in pseudospin 
space. 

As is familiar from the case of quantum Hall ferromag- 
nets^S such a collective state is represented as a texture 
state 



\G) =e-*°\G)\ n , 
where the operator e~ with 

° = E^p lf7 P 5 -. 



(3.10) 



(3.11) 



locally tilts the pseudospin from n to n p by small an- 
gle n p ~ n x n p . Repeated use of the charge alge- 
bra (|2.13[) then allows one to express the texture-state en- 
ergy (G\H\G) = (G\e i0 He- i0 \G) in terms of the struc- 
ture factors in Eq. (|3.4p . and this yields an effective 
Hamiltonian as a functional of Q p or its x-space rep- 
resentative fi a (x, t). 

The angle variables fi a (x, t) may be normalized so that 
^ a (£! a ) 2 = 1 classically. As we shall see below, the ef- 
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fective theory is expressed in terms of the following com- 
ponents of O a , 



V 

C 



O 1 cos 9 -O 3 smt 

o 2 , 



(3.12) 



along the two orthogonal axes (i.e., the tilted " 1" axis and 
the "2" axis) perpendicular to the pseudospin n, with 
normalization (n a U, a ) 2 + rf + £ 2 = 1. (Here we have 
chosen = and set {n a } = (s'm9, 0, cos#), as remarked 
above; the </> ^ case, if needed, is readily recovered^) 
Actually they refer to the following induced pseudospin 
components in the excited state |G), 

r, ~ -^- 1 (G\S 2 \Cf), 

C ~ 7 p 1 (G|S 1 cos6>-S 3 sin(9|G), (3.13) 

as seen from the induced pseudospin 6(G\Sp\G) » 
( Po /2) lp e abc n b p n c . 

We expand (G\H\G) to second order in O and retain 
all powers of derivatives acting on O to study the spec- 
trum over a wide range of wavelengths. The calculation 
is outlined in Appendix B. The result is 



(G\H\G) = 


p Q J d 2 x(H cS +Wcoll), 




Wool] = 




- n ch + sn, 


r, = 


-£/sin6» + ^ p , 




r< = 


E"(6) + G p , 




sn = 


£'(0){C + OV(2sin0)}; 


(3.14) 



E'{6) = dE{9)/d9, etc. Here 



G r 



Wp = 



_£ 

2 

_£ 
2 

2 



1 TT 
2V2 

1 /tt 

2V2 



„ PxPy . „ £p x 

2 — tt A„ cos 9 + i — — t sm f 



6? + P- \ 



£, q — b q sin 2 9 — P- X q cos 2 



V2' 9 

with V/ c = a/(e h £), P_ = (p 2 -p 2 )/p 2 and 

/» OC 1 \ 

C = ^^ /2 f-/ ^^ 2/2 (i-^ 2 ) JoN) 



(3.15) 



2 2. /0 



dze z 1 2 z 2 J 2 (zq), 



b q = ve-i 2 / 2 (q/2-q 3 /8) 

+ - / dze- z ^ 2 (l-4z 2 + z A )J (zq), 

r = V e~^' 2 q/2 
2 
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2 / 2 z(l-\z 2 )j'{zq), 



(3.16) 



where q = £|p|; substitution p — ► — zV is understood in 
the x representation. For p — > 0, £ — > —(1/2) y/n/2, 



(X q , b q ) — > and r 9 — > —(1/4) y/ir/2 while they all tend 
to zero for p — > 00. See Appendix C for more explicit 
forms of the integrals involved in these functions. 

In Hcoii, Tich refers to a topological charge (to be de- 
tected with a constant potential ^4q) 



3 Jd 2 TiH c h = -e Jd 



r, l^xh,u= -< /.PxAi-e afc e < j(3 i fl )(9 j n i ) ) i c , 

(3.17) 

where e xy — 1 and e 123 = 1. With n c promoted to O c , 
J <i 2 x7i c h involves the winding number, which implies 30 
that possible topologically-nontrivial semiclassical exci- 
tations (Skyrmions) associated with O a , in general, carry 
electric charge of integral multiples of v e. Note also that 
Sn, involving a term linear in £, disappears for a stable- 
state configuration for which E'(9) = 0. 

The the kinetic term for O a is supplied from the elec- 
tron kinetic term as Berry's phase, 



(G\id t \G) = -^p j2 n ^ abCQ -k d t n i 

k 

Jd 2 x Cd tV , (3.18) 



Pn 
2 



to 0(00), apart from surface terms. This shows that £ 
is canonically conjugate to (po/2)rj. One can now write 
the effective Lagrangian for the collective excitations as 



L = 2 C <9t ?7- Wcoll [?/,C]- 

Note here that this Lagrangian is written as 
po [d 2 *L= (G\e w (id t - H)e~ l °\G). 



(3.19) 



(3.20) 



This representation realizes and systematizes the single- 
mode approximation 29 (SMA) within a variational 
framework* 3 ^ The present theory thus embodies nonper- 
turbative aspects of the SMA. 

Upon elimination of (, Eq. (|3.19[) leads to an alterna- 
tive form of the effective Lagrangian for r\ as follows: 



1 



1 



U = -(W--*(M P )'$, 



m p = 2 x /r,r c -|H/p| 2 , 



(3.21) 



where we have set $ = (1/2) (r^) -1 / 2 r\. 

The spectrum M p of collective excitations is in general 
anisotropic in p at low energies and depends critically 
on the stable-state configuration n. In particular, the 
leading long-wavelength correction in n co \\ starts with 
the direct Coulomb interaction of O(p), which leads to 
the spectrum, 



M, 



J{2K n K C ) 2 +Vc p{« 2 COS 2 9 p 2 y + K 2 P 2 x } 7 

r^| P =o = -£/ sin6>, 



Clp=o 



E"(9). 



(3.22) 
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The excitation gap at zero wave vector is thus given by 
M p _»o = 2k i? k^. In contrast, M p recovers isotropy and 
the standard excitonic behavior— at short wavelengths, 

^oo«V r i/2 + /t*+4 (3.23) 

with T,, -> Vi/4 + «2, r c -> Vi/4 + k 2 c and W p -> for 
p_>oo; Vi = V cy /V/2. 

It is worth noting here that the Coulomb interaction 
alone yields k 2 — 0, i.e., a flat direction in (77, Q space. 
This implies that, unlike in ordinary bilayer QH systems, 
there is no cost of interlayer capacitance energy for the 
pseudo-zero-mode sector which essentially resides in the 
same layer. Coherence is thus easier to form in this sector 
of bilayer graphene. 

If we set — » 0, our M p precisely reproduces the exci- 
tation spectrum derived in Ref. [26j by assuming spatial 
isotropy; actually, the effective Hamiltonian of Ref. [26| 
is apparently different from our 7Y co ii but the spectrum 
turns out to be the same. It is our use of general tex- 
tures n that allows 7i co \\ to handle spatially anisotropic 
situations as well. 



IV. POSSIBLE GROUND STATES AND 
COLLECTIVE EXCITATIONS OVER THEM 

In this section we study the spectrum of the half-filled 
pseudo-zero-mode state and the associated collective ex- 
citations. Let us first gain a rough idea of the strengths 
of to and E|| relative to the Coulomb correlation energy 
AE C . A naive estimate 



A£ c = -^-- Rj 2.2 v ^[T]meV (4.1) 

with a typical value e& ~ 4 indicates that the bare 
Coulomb interaction is apparently sizable, as compared 
with the basic Landau gap lo c ~ 3.9 B[T] meV. We re- 
mark that Eq. (|4.ip is likely to overestimate AE C . Actu- 
ally AE C is written as an integral of the form 

AE C ~ 2]T« p e-« 2 / 2 ( g 2 /4)(l - q 2 /A) (4.2) 
p 

with q = £\p\, and the main contribution comes from the 
momentum region \p\£ ~ 1, where the Coulomb interac- 
tion is very efficiently weakened [v p — > v p /e(p)], as indi- 
cated by a random-phase-approximation study^ of the 
dielectric function e(p). This screening effect essentially 
comes from vacuum (Dirac-sea) polarization, specific to 
graphene. It may effectively be taken care of by setting 
(b — * £fc£sc; a simple estimate™ gives e sc ~ 9 at B ~1T 
and e sc - 3.6 at B ~10T. 

The ratio of the intrinsic zero-mode level gap 2m = 
0.013 B[T] U to AE C is generally small, 

m/AE c w 3e sc x 10~ 3 y/B[T] E/[meV], (4.3) 
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FIG. 1: Excitation spectra M p , plotted in units of V c = a/e^l 
for m/AE c = 0.01, 0.05, 0.13, 1.1 and 1.2. The real curves 
refer to the profiles in p x at p y =0 (i.e., normal to the in-plane 
field E y ) and dashed curves to those in p y at p x = 0. 

for a band gap U of 0(1 meV), and increases with B and 
U. For the in-plane field £ — e£\E^\/^/2 the ratio 

£/AE c « 0.9 e sc x 10~ 3 E[V/cm\/B[T] (4.4) 

is on the order of 10% for E = |E|| | = 10 V/cm at B = IT 
with e sc ~ 9. Note that £ « to for £ w 3.6 V/cm x 
U[meV] B[T] 3 / 2 . 

To determine the orbital configuration of the half-filled 
zero-mode state \G) one has to look for the minimum of 
E{6), E min = E{8 min ) with E'(6 min ) = 0. 

Let us begin with the case where Ey is absent. It is 
clear from E(9) of Eq. (|3.9p that the Coulomb correlation 
favors 9 = while the intrinsic asymmetry m > alone 
favors 9 = tt; AE c and m thus compete. 

(i) In case m > AE C (although rather unrealistic), one 
finds E m i n = —(2m — AE C ) and k 2 = k 2 = m — AE C at 
9 = tt. The collective excitations have a finite energy gap 
M p=0 = 2 (m — AE C ) > and the spectrum is isotropic, 
reflecting the rotational invariance of the bilayer system 
about the z axis || B, or more explicitly, invariance of 
AH + H C (for A = 0) under 17(1) rotations in Eq. (pH]) . 

(ii) On the other hand, for < m < AE C , one 
finds E min = -AE c (m/AE c ) 2 at 8 = ±6> min with 
sin (#min/2) = m/AE c ; 9 m i n varies from to tt with 
increasing m. Here we encounter a somewhat strange 
situation. Note that, for 0^0, ±7T, the rotational invari- 
ance is spontaneously broken. Indeed, one finds k 2 = 
and k 2 - = 2m (1 — m/AE c ). Accordingly the collective 
excitations about this state are gapless (as the Nambu- 
Goldstone modes) , and the spectrum p.22j) is anisotropic 
in p. 

The excitation spectrum M p in general exhibits a lo- 
cal minimum (roton minimum) around \p\t ~ 2; see 
Fig. 1 (a). In the < m < AE C range of case (ii), the ro- 
ton minimum changes critically with to: The minimum in 
p x comes down, as m is increased from zero, and touches 
zero (gap) at m/AE c w 0.131 (with 6> min w 42.6°). The 
spectrum loses sense (becoming pure imaginary) until m 
reaches m/AE c w 0.869 (with 6» min w 137.4°) where the 
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roton minimum reappears; it then rises and returns to the 
the m = spectrum at m — AE C . This peculiar feature 
inspires one to find an interesting structure of E{0), 

E(9;m) = E(tt - 9; AE C - m), (4.5) 

valid for £ ^ as well. This implies that the half-filled 
state realized at 6* m in for a given m < AE c /2 is paired 
with the state realized with angle 7r — 9 m i n at a larger 
intrinsic gap oc m! = AE C — m. Both the energy E(9 m i n ) 
and excitation spectrum M p completely agree for this 
pair of states, as seen from Eqs. (|3.15|) and (|3.21[) . In 
particular, the filled n — level realized at m = and 
the filled n = 1 level realized at m = AE C share the same 
energy and collective excitations. We thus find a kind of 
(small-m/large-m) duality in the < m < AE C range. 

In this m range the texture state, taken to be homo- 
geneous in space, acquires spontaneous in-plane electric 
polarization oc sin# m i n . The anomalous behavior of the 
roton spectrum about m/AE c ~ 0.5 or 9 ~ 7r/2, men- 
tioned above, reflects a potential instability of the tex- 
ture state due to spontaneous polarization. Such electri- 
cally polarized homogeneous configurations, unless po- 
larization is relatively weak, are unstable against local 
charge inhomogeneities and would decay into inhomoge- 
neous configurations. 

A local charge excess would align electric dipoles out- 
ward or inward, and let them drift in a magnetic field. 
One may thus imagine a picture of charged electric 
dipoles drifting around local charge centers (distributed 
randomly or in some patterns on the real sample and 
substrate), clockwise or anticlockwise depending on the 
sign of the local excess charge. We speculate that the 
half-filled state in the realistic < m < AE C range may 
form many such domains for stabilization. 

Let us next set m — > 0, i.e., consider the case of zero 
band gap U — and study the effect of E||. The in- 
plane field £ = e£Ey/y/2 > tilts the pseudospin toward 
9 = —n/2 and competes with AE C which favors 9 = 0. 
As a result, 9 m { n varies from to —n/2 as E y is increased. 
The charge carriers thereby acquire a nonzero electric 
dipole moment oc sin 9 and the pseudospin waves always 
have a finite excitation gap for £ ^ 0; see Fig. 2. For weak 
field A£/AE C = R < 1, one finds 6> min « -i? 1 / 3 , E min « 
-(3/16)i? 4 / 3 A£ c , and n\ « Zk% m (3/4) i? 2 / 3 AE C so 
that the excitation gap grows as 

M p=0 m (V3/2) (A£/AE C ) 2/3 AE C . (4.6) 

For larger £ ^> AE C the gap rises almost linearly with £, 

M p=0 ^2£ + AE c /2, (4.7) 

along with k 2 w £ and k 2 w £ + AE c /2. As seen from 
Fig. 2, Eqs. (|4.6[) and (|4.T[) combine to give a practically 
good description of the excitation gap over the entire 
range of £ ; crossover takes place around £ / AE C ~ 0.3. 
Note that the Coulomb correction significantly enhances 
the excitation gap; in particular, the gap rises promi- 
nently as £ 2 / 3 for a weak field. 
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FIG. 2: (a) Excitation gap at zero momentum, plotted in 
units of AE C as a function of the in-plane field, £/AE c . The 
red dashed curve refers to Eq. (|4.6[) and the green dotted line 
to the asymptotic form of Eq. (|4.7|l . (b) Angle of inclination 
#min in degrees, (c) Excitation gap M p= o, in units of meV and 
GHz, plotted as a function of the in-plane field \E« \ for some 
typical values of (S[T], e sc , [/[meV]). Dashed curves refer to 
the cases where the screening effect is turned off, esc — * 1. 

For £ 7^ the excitation spectrum M p is necessarily 
anisotropic, especially at low momenta, and, as seen from 
a contour plot of M p in Fig. 3, anisotropy of the roton 
minima already develops around £/AE c ~ 0.02. The 
spectrum recovers isotropy for \p\£ ~ 5 or larger and the 
asymptotic spectrum is lifted roughly by the amount of 
the excitation gap, as seen from the spectrum profiles in 
Fig. 3(a) and from Eq. (|3.23|) . Note that the roton mini- 
mum, unlike the m ^ case, shows no sign of instability. 

When both m and £ are present, their effects generally 
tend to add up. The texture excitations always have 
a gap and their potential instability, in a certain range 
about m/ AE C ~ 0.5, weakens and eventually disappears 
with increasing £; the roton dip remains to be a local 
minimum as long as £ > m. We remark that the duality 
implied by Eq. (|4.5|) would hold in the presence of £ as 
well. 

In Fig. 2(c) we plot the excitation gap as a function of 
|E||| for some typical values of B[T] and m oc [7[meV]. 
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FIG. 3: (a) Excitation spectra M p , in units of V c = a/tbl, 
for £/AE c = 0.01, 0.05 and 0.1. The real curves refer to the 
profiles in p x and dashed curves to those in p y . (b) Contour 
plot of M p for £/AE c = 0.02. 



With the effect of screening e sc taken into account, M p= q 
falls in roughly the same frequency range of microwaves. 

If E|| is sufficiently strong, a sizable excitation gap 
M p= o would arise, leading to an incompressible v = 2 
state. One could thereby observe the (spin-degenerate) 
v = 2 Hall plateau with a suitably strong injected current 
for U ^ and U — as well. 

Of the pseudospin cx n the angle 8 is related to the 
ratio in amplitude of the n — 0+ and n = 1 modes while 
the angle <j) is related to the relative phase between them. 
One would now have control of mixing of the zero-mode 
levels by adiabatically changing the strength and direc- 
tion of an injected current. 



V. CURRENT AND RESPONSE 

In this section we study the response of the pseudospin 
waves in the (0+, 1) sector under uniform fields B and 



E y to a weak time- varying external field. To this end we 
consider a weak vector potential A(t) = (A x ,A y ), which 
describes an external field ~ dt A and, at the same time, 
serves to probe the current. For simplicity we take it to 
be spatially uniform. 

The current operator SH/5A derives from two parts 
in H of Eq. (|2.1[) . One coming from Ho is the ordinary 
form of current, which has no projection to the (0+, 1) 
sector and induces transitions to other levels. The other 
one, specific to bilayer graphcnc, derives from the O(zU) 
part of H as , and the relevant O(A) portion of its (0+, 1) 
projection is written as 

H A = -2V2mp e£{A x Sl =o -A y Sl =o ), (5.1) 

where m = zU/A. Evaluating {G\Ha\&) yields an addi- 
tion to H co \\ of the form 

Ha = -V2e£{coA x +c x A x ( + c y Ayi]}, 
co = msmO, c x = mcos(9, c y = m. (5-2) 

One might read from Ha the current density carried 
by the zero-mode sector as j x m = —y/^Poeim sin^min 
[which, for V c = 0, equals — (ye 1 jh) E y ]. This, unfortu- 
nately, is not a complete amount of current yet. Inter- 
Landau-level transitions caused by the ordinary current 
induce some extra charge in the (0+, 1) sector. In other 
words, the presence of A causes level mixing, which mod- 
ifies the current within the (0+, 1) sector. Such a modifi- 
cation was calculated earlier— for standard QH systems 
and, in the present case, it is given by 

5p- P « (S°- P -S 3 _ p ) u u (p), 
Ull (p) w 2ief (1 - £ 2 p 2 /4)( Px A y - Py A x ), (5.3) 

to O(A), apart from terms of 0(d t Ai/id c ). This itn(p) 
represents charge accumulated in the n = 1 level via the 
n = 1 — * ±2 — > 1 intcr-Landau-level transitions. 

The induced charge also carries current within the 
(0+, 1) sector via the interaction 

5 Ha = ^ v p : {ft" S P-v} : ~ e ( A o) P § P-p- ( 5 - 4 ) 
p p 

The current response is now calculated from (G\Ha + 
SHa\G). The result again takes the form of Eq. (|5.2[) . 
with coefficients (co,c x ,c y ) modified as follows 

c = m sin 9 + (1 - cos 6) £ + (Vi/32) sin 26, 

= E-E'(6), (5.5) 

c x = mcos6 + --- = -E"(9), (5.6) 

C y = m + (Vi/16)(cos0- 1), 

= £ cos / sui6 - E' [6)/ sinO. (5.7) 

With E'(9) — > 0, Eq. (|5.5|) verifies that the pseudo-zero- 
mode sector carries the correct amount of Hall current (in 
response to a uniform field £ ) with conductance a xy = 
-ve 2 /h. 
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The Hamiltonian Ha + SHa governs the microwave 
response of the pseudo-zero-mode sector. It is combined 
with H C oii to yield the source term <I> {2c y (T^Y^Ay — 
M r c)~ 1/2 d t A x for L$ in Eq. (|33T|) . Solving for the 
stationary action then yields a response of the form ~ 
A y {- ■ -)dtA x . From this one can read off the optical Hall 
conductance due to virtual transitions within the (0+, 1) 
sector, 



v e 



Aa xy (u;) — - — — cos 9 
ZttTl 



Ml 



Ml 



(5.8) 



which is significantly peaked around lu ~ M = M p= o, 
the pseudospin-wave gap. 

The collective excitations within the (0, 1) sector thus 
contribute the cos 9 portion [Aa xy (u) — > 0)] of the Hall 
conductance a xy while the remaining (1— cos 9) portion of 
<j xy essentially comes from the n = 1 — * ±2 — * 1 virtual 
transitions with larger gaps ~ \[2lo c ^ Mo. These two 
components are distinguishable via microwave or light 
response. 

With disorder taken into account, the diagonal con- 
ductivity Aa xx (uj) also is significantly peaked around 
u) ~ Mq 1 and Mq varies critically with £ or by an in- 
jected current, as we have seen in Fig. 2. Microwave or 
infrared experiments,— via absorption, reflection or con- 
ductance fluctuation, would provide a direct means to 
explore such unique dynamics of the pseudo-zero-modes. 



VI. SUMMARY AND DISCUSSION 



zero-modes (i.e., relative phase and magnitude) is also 
externally controllable to some extent. 

The pseudo-zero-mode sector of bilayer graphene is 
particularly suited for exploring coherence phenomena. 
This is because it essentially resides on the same layer so 
that, unlike in ordinary bilayer QH systems, there is no 
cost of inter layer capacitance energy for it. 

An experimental signature of the field-induced gap is 
to observe the quantum Hall effect with an injected cur- 
rent; one would be able to resolve the v = ±2 Hall 
plateaus (or the spin-resolved v = ±1 plateaus) using 
a suitably strong current. 

The collective excitations within the pseudo-zero-mode 
sector also carry a considerable portion of the total cur- 
rent. A direct study of the excitation gap A/ p= o an d 
its field dependence by microwave absorption or reflec- 
tion would clarify the unique controllable features of the 
pseudo-zero-mode sector in bilayer graphene, and, in ad- 
dition, the effect of screening on the Coulomb correlation 
energy AE C due to vacuum polarization. 
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Zero-mode Landau levels, specific to graphene in a 
magnetic field, are very special. Their presence has a 
topological origin in the chiral anomaly. They show quite 
unusual dielectric response that reflects quantum fluc- 
tuations of the vacuum state (the Dirac sea). Bilayer 
graphene supports eight such zero-mode levels which, 
unlike in monolayer graphene, involve two different or- 
bital indices n = 0, 1. As a tunable band gap develops, 
four of them at one valley are isolated from the others 
at another valley and remain nearly degenerate although 
its fine structure sensitively depends, via mixing of zero- 
modes, on the environment. 

In this paper we have studied the effects of an exter- 
nal field and the Coulomb interaction on such an isolated 
zero-mode quartet. This pseudo-zero-mode sector, espe- 
cially at half filling, supports, via orbital mixing, quasi- 
particles with charge and electric dipole, which give rise 
to characteristic collective excitations, pseudospin waves. 
We have constructed a low-energy effective theory of 
pseudospin waves with general pseudospin textures and 
noted a duality [Eq. (|4.5[) ] in the excitation spectrum. 
The excitation gap at zero momentum turns out to be 
generally small, reflecting the intrinsic degeneracy of the 
pseudo-zero-mode sector, and the Coulomb exchange en- 
ergy works to enhance the effect of the in-plane field on 
the gap. This means that the gap is tunable by an in- 
plane field or by an injected current; the mixing of the 



APPENDIX A: STATIC STRUCTURE FACTORS 

In this appendix we outline the derivation of the pseu- 
dospin static structure factors (<Sp5 q ) [in Eq. (|3.4[) ] for 
the half- filled pseudo-zero- mode levels \G) with (Sq) = 
(N e /2) n a pointing in a general direction n = {n a } in 
pseudospin space. Let us first note that, when only 
the n = 0+ level is filled, i.e., for n 3 = 1 polariza- 
tion, these structure factors are readily calculated: Fill- 
ing the n = 0+ level and leaving the n = 1 level empty 
immediately imply that (SgS^) — (po/2) 2 <5 p ,o <^q,o for 
fi,u e (0,3), and (: S£S q :) = for a, b e (1,2). One 
may then note the algebraic relation 



+ ^e abc S c p+^ ) + 



°p°q 



(Al) 



with / pq = i 2 (p ■ q — ip x q)/2, and determine, e.g., 
(: S p S 3 _ p :> = (iY e /4)( Po <W - 7*) and (S p Sl p ) = 
(A e /4) 7 2, with 7p = e-'V/4. 

The structure factors for the half-filled state with 
a general pseudospin polarization n are obtained from 
these n 3 = 1 structure factors by a suitable rotation in 
pseudospin space. Note first that the n 3 = a 3 = ±1 
eigenspinors | ± 1) of the Pauli matrix cr 3 are rotated by 
angle (9, (f>) to form the <j a = ±n a eigenspinors U\ ± 1) 
with U = e' 1 ^" 7 I" 2 e~ l6a I" 1 . Accordingly we decompose 



10 



the zero-mode field <f = (ipo+^iY [defined in Eq. I|2.6[) ] 
into the a a = ±n a eigenmodcs by writing \P = U^' . 

On substitution \& = [/'I'', Sp 1 are rewritten as lin- 
ear combinations of the pseudospin operators S' p a ~ 
*' t (CT a /2)e 4p r *' composed of S° = Sg. One can 
then calculate (S£S^) for general n a from the structure 
factors (S'^Sq) for the n 3 = 1 state. The result is sum- 
marized in Eq. (|3.4|) . 

Note, in particular, that S£ = n a S® + ■■■■ This 
yields = n°(5^ 3 ) = n a (p /2) <5 Pi0 and tells us that 
the normal-ordered factors take particularly simple form 
(: =n«n 6 <:S*S*:>, or 

(:S^ p :> =-n^(iV e /4)( 7 2-p 5 p>o ), (A2) 

with n° = 1, as quoted in Eq. (|3.6p . 



^2 = -^^ra^/pe^e^s 1 ^ 1 , 

iv 3 = n a _ k n^\w° p \ 2 s ^ a \ 

iV 6 = n a _ k n b k e abc w p w^ p s^ cl3 \ (B3) 
where k x p = k x p y — k y p x . 

One can evaluate J2 P v pJp by integrating over p and 
leaving the k integration as it is. One may express the 
sines and cosines in terms of e ±lt pxk . Integration over 
p is then carried out as a Fourier transform of the form 
^ p u p e ~ l p / 2 +*p x (powers of pi) with x — > £ 2 (xk) = 

^ 2 (ky > 

The rest of terms in 7i co ii are obtained via the induced 
pseudospin to 0(£lil) 



APPENDIX B: COLLECTIVE EXCITATIONS 

In this appendix we outline the derivation of TL co \\ 
in Eq. (13. 14|) . Let us first consider the contribution 
from the Coulomb interaction, | Tl p v pJp with J p = 
(G\(p p )° (p- P )°\G) and (p p )° = e i0 p p e- i0 . Expand- 
ing (p p )° in powers of O oc f) by repeated use of the VFoo 
algebra ()2.13|) and subsequently substituting the struc- 
ture factors in Eq. (|3.4|) allow one to evaluate J p . 

The 0(17) term is thereby written as 



J« = -porin^oe^w 



po%%{\K\ -K\ )* 



,3|2^ .{1,3} 



(Bl) 



a c/3 . 



under a symmetric integration over p; s 
for short. (Here, we have employed the convention 
and n 2 = 0, as remarked in the text.) 
Similarly, the 0(£!£1) term is written as 



7(2) 
J P 



~2 _,2 |_a5c„a„„6 oc |2 

Po7 P |e " W p"- P l 

+f 7^E cos2 (^ 2kx p) 



-y7pT sin2 (^ 2kx p) (^ + ^4) 

k 

-^ 7p ^sin(^ 2 kxp)(2iV 5 + 7V 6 ), (B2) 

k 



with 



b OiB-y abc A^i,c\ 



a\0\ 



Po 
2 



n a 5 P , + e abc % n 



{n a (n b ,n b ) p -(n a ,n») p n b } , (B4) 



where (fi a , n b ) p = £ k cos(£ 2 kxp/2) k+p « k for short. 
This, in particular, is used to evaluate the contribution 
(A°) from A in Eq. (f2~T2l . Somewhat tedious cal- 
culations along these lines eventually lead to 7i C oii in 
Eq. (I3~T41 . 



APPENDIX C: INTEGRALS 



The integrals appearing in Eq. (|3.16p , apart from their 
overall factors, are expressed in terms of the modified 
Bessel functions 



dze-^l 2 (---)=e-t 2 l A J^ [c I (q 2 /A) + c 1 I 1 (q 2 /4) 



with coefficients 

(co,ci) = 

(co,ci) = 

(co,ci) = 

(co,ci) = 



(CI) 



(l/4)(2 + g 2 ,-g 2 ) for£„ 
(l/2)(q 2 ,-2-q 2 ) for A„ 
(g 2 /T6)(-2 + <? 2 ,-<? 2 ) for b q 
-(q/8)(l + q 2 ,-3-q 2 



for r q . (C2) 
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